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ABSTRACT 

The galaxy CXOC J100043. 1+020637, also known as CID-42, is a highly unusual object. 
An apparent galaxy merger remnant, it displays signatures of both an inspiraling, kiloparsec- 
scale active galactic nucleus (AGN) pair and of a recoiling AGN with a kick velocity > 1300 
km s~^. Among recoiling AGN candidates, CID-42 alone has both spatial offsets (in optical 
and X-ray bands) and spectroscopic offsets. In order to constrain the relative likelihood of 
both scenarios, we develop models using hydrodynamic galaxy merger simulations coupled 
with radiative transfer calculations. Our gas-rich, major merger models are generally well 
matched to the galactic morphology and to the inferred stellar mass and star formation rate. 
We show that a recoiling supermassive black hole (SMBH) in CID-42 should be observable 
as an AGN at the time of observation. However, in order for the recoiling AGN to produce 
narrow-line emission, it must be observed shortly after the kick while it still inhabits a dense 
gaseous region, implying a large total kick velocity (uk ^ 2000 km s^^). For the dual AGN 
scenario, an unusually large broad-line offset is required, and the best match to the observed 
morphology requires a galaxy that is less luminous than CID-42. Further, the lack of X-ray 
emission from one of the two optical nuclei is not easily attributed to an intrinsically quiescent 
SMBH or to a Compton-thick galactic environment. While the current data do not allow either 
the recoiling or the dual AGN scenario for CID-42 to be excluded, our models highlight the 
most relevant parameters for distinguishing these possibilities with future observations. In 
particular, high-quality, spatially-resolved spectra that can pinpoint the origin of the broad 
and narrow line features will be critical for determining the nature of this unique source. 
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1 INTRODUCTION 

Recent breakthroughs in numerical relativity have allowed for fully 
general-relativistic simulations of black hole (BH) mergers, in- 
cluding calculations of gravitational waveforms and kick velocities 
from gravitational-wave (GW) recoil (e.g., |Pretorius|[2005| |Cam-| 
Ipanelli et al.|2007a| |Baker et al.|2008[ l. Surprisingly, these simu- 
lations have foun d that GW recoil kicks may be quite large, up to 



5000 km s"^ ( |Campanelli et al.||2007b 



Lousto & Zlochower 



[20TT] >. Thus, the merged supermassive black hole (SMBH) may be 
ejected from its host galaxy in some cases, and kick velocities of 
even a few percent of this maximum could be of astrophysical rel- 
evance. This discovery has prompted many studies of the possible 
electromagnetic (EM) signatures from recoiling SMBHs as well as 
observational searches for such objects. 

While prompt EM signatures of SMBH mergers are of great 
interest for identifying counterparts to eventual GW detections, 
long-lived signatures such as offset AGN provide the best chance 
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of finding recoiling SMBHs using purely EM observations. An 
AGN ejected from the center of a galaxy can carry along its ac- 
cretion disk as well as the broad-line region, which could produce 
an AGN spatially offset from its host galaxy or an AGN with offset 
broad emission lines (Madau & Quataert 2004; Loeb|2007[|Blecha| 
[& Loeb 2008; Komossa & Merritt 2008). In either case the offset 
AGN lifetimes may be up to tens of Myr for a fairly wide range in 
kick speeds ( [Blecha et al.|201 1) . To date, several candidate recoil- 
ing SMBHs have been found with either spatial offs ets jBatcheldor 



|et al. 2010, Jonke r et al.|20 10|l or kinematic offsets jKomossa et al. 



2008, Shields et al. 2009 , Robin son et al.|2010| l, though none have 
been confirmed ( Komossa,2012. for a review, see). 

Recent years have also seen renewed energy in the search 
for inspiraling pairs of SMBHs. Systematic searches revealed that 
about 1% of AGN have double-peaked narrow [OIII] lines, which 
could be a signature of SMBH orbital motion on kiloparsec scales 
( IComerford e"raL]|2009a| |Liu et aI.||20T0bl [Smith et al.|[20T0l ). 
Follow-up observations using higher-resolution imaging and slit or 
integral field unit (IFTJ) spectra have identified many strong dual 
AGN candidates ( |Liu et al.|[20T0al |Shen et al.|[20ria| |Fu et al.| 
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|2011a| [20121 IMcGurk et al.|2011|[Rosario et al.|20lTl|Comerfordl 
|et al.||2012[l, and in tw o cases, evidence for dual compact X-ray 
( [Comerford et al.|201 lb and radio ( Fu et al.|201 lb| l sources, respec- 
tively, has been found. Shen et al. ]( |20T lab estimate that at least 10% 
of double-peaked NL AGN may actually contain dual SMBHs. 

Among these new data and interpretations of SMBHs in merg- 
ing galaxies, the galaxy CXOC J100043. 1-1-020637 (also known as 
CID-42; z = 0.359) is a fascinating case. Discovered serendipi- 
tously by [Comerford et aL] ( |2009bl hereafter JC09) based on the 
striking morphology evident in its COSMOS HST/ ACS image, it 
is an apparent face-on spiral with a prominent tidal tail and two 
well-resolved brightness peaks in its center, separated by 2.46 kpc 
(Figure [T^). The disturbed morphology strongly indicates a recent 
galaxy merger. DEIMOS slit spectroscopy revealed dual emission 
components in three different narrow lines (NLs) with a velocity 
separation of ~ 150 km s^^ and a spatial offset roughly consis- 
tent with the separation of the two brightness peaks in the image, 
albeit with low signal-to-noise (S/N, JC09). No evidence for broad 
lines (BLs) is apparent in this spectrum. Thus, JC09 conclude that 
CID-42 is a good candidate for a kpc-scale. Type 2 AGN pair. 

Around the same time, CID-42 was also identified as an un- 
usual object during visual inspection of the optical counterparts to 
Chandra-COSMOS sources ( [Elvis et al.|2009l[Civano et al.|2012a| l. 
Additionally, a broad (H/3) line was detected in the SDSS spectrum 
of CID-42 fShen et al. 201 lb I that was not seen in the DEIMOS 
spectrum. jCivano et al. (2010, hereafter EC 10) and [Civano et al.[ 
( [2012b| hereafter EC 12) also found a prominent broad H/3 line 
in IMACS/Magellan CTrump e t al.|2007j[2009) l and VIMOSA^LT 
(2COSMOS; [Lilly et al. 2007]l spectra. More interestingly, the 
broad and narrow H/3 components are ojfset by ^ 1360 km s~^. 
The double-peaked narrow lines seen in the DEIMOS spectrum 
would be unresolved in these spectra, which have higher S/N but 
lower spectral resolution (R=700 on IMACS and VIMOS versus 
R=3000 on DEIMOS). There is also no measurable spatial offset 
between the broad and narrow emission components, though this 
may be attributable to seeing conditions (FCIO). 

A BL/NL offset of 1360 km s^^ is much too large to be ex- 
plained by orbital motion of dual AGN on kiloparsec scales. How- 
ever, this is a signature expected from a recoiling AGN, which 
should retain only the most tightly bound gas and stars (including 
the BL region) and leave behind gas and stars at larger radii (includ- 
ing the NL region) ( [Merritt et al.|2006y Accordingly, EC 10 pro- 
posed that CID-42 may contain a single, rapidly-recoiling SMBH 
rather than two inspiraling SMBHs. In this scenario, one of the op- 
tical brightness peaks (the NW source in the HST image. Figure 
[T^) corresponds to the central stellar cusp left behind by the recoil- 
ing SMBH, and the other (SE source) corresponds to the accretion 
luminosity of the recoiling AGN. A surface brightness decomposi- 
tion of this image performed by FCIO is consistent with this pic- 
ture, in that the SE source (the purported recoiling AGN) is best fit 
as a point source, while the NW source (the purported empty stellar 
cusp) is best fit with an extended Sersic component. 

Even more intriguing evidence as to the nature of CID-42 
comes from High Resolution Chandra observations (FC12), which 
reveal that the X-ray emission detected in CID-42 arises from only 
the SE source, with a 3(j upper limit of 4% of the observed X-ray 
flux from the NW source. This finding supports the recoil scenario 
for CID-42, in which only one of the optical brightness peaks con- 
tains an AGN, but it does not rule out the possibility of a quiescent 
or Compton-thick AGN in the NW source. 

Another possibility is a gravitational slingshot recoil result- 
ing from a triple-SMBH interaction. This may occur if the binary 



SMBH inspiral time is long enough that the host galaxy undergoes 
a subsequent merger in the meantime, bringing a third SMBH into 
the system. We focus instead on GW recoil in our modeling, on 
the grounds that a triple-SMBH interaction is most likely in a gas- 
poor merger remnant, unlike CID-42. Nonetheless, a slingshot re- 
coil cannot be ruled out by the observations, and this possibility is 
discussed further in § |3.1| 

Because CID-42 displays signatures of both dual and recoil- 
ing AGN, and because it is the only recoil candidate with evidence 
for both spatial and kinematic offsets, this object is an ideal case 
study for theoretical models. [Blecha et ar] ( [201l | have conducted a 
large parameter study of recoiling SMBHs in hydrodynamic simu- 
lations of galaxy mergers, which allowed them to estimate observ- 
able lifetimes of spatially- and kinematically-offset recoiling AGN 
for a wide range of kick velocities and merging galaxy parameters. 
Here, we use hydrodynamic simulations to constrain the possible 
parameters of a recoiling or dual AGN in CID-42 and to gain in- 
sight about the relative plausibility of these scenarios. We provide 
models for both a GW-recoiling AGN and for an AGN pair, based 
on the methods of [Blecha et"aL[ ( |2011^ and [Blecha et al.[ ( |2012| l. 

The details of our simulations are outlined in § [2.1[ We de- 
scribe the properties of our best-fit recoiling SMBH simulations in 
§ [2.2[ and those of our best-fit dual SMBH simulation in § [2.3[ In 
§ [3.H3.3[ we discuss the plausibility of the recoiling and dual AGN 
scenarios, and we briefly summarize in § [3.4[ 



2 MODELS 

2.1 Initial Conditions and Merger Simulations 

To model the possible origins of CID-42, we begin by simulating a 
major galaxy merger using the smoothed-particle-hydrodynamics 
(SPH) code GADGET-3 (Springel 20051. Each progenitor galaxy 
contains dark matter, star, and gas particles, as well as an accreting, 
central SMBH. GADGET-3 also contains sub-resolution models for 
star formation and for supernova and AGN feedback. 

We have designed simulations for both the GW recoil and dual 
SMBH scenarios that match as closely as possible the observed 
properties of CID-42. Specifically, the progenitor galaxies have a 
mass ratio of g = 0.5, and the primary galaxy has a baryonic (total) 
mass of 1.9x 10^^° Mq (4.5 x lO'^^ Mq). Because star formation oc- 
curs throughout the simulation, the initial baryonic mass is chosen 
such that by the time the system would be observed, near the time 
of SMBH merger, the stellar mass matches that inferred for CID-42 
(~ 2.5 X 10^" M0, FCIO). The galaxies are disk-dominated with 
10% of the initial stellar mass in a bulge component. For the recoil- 
ing (dual) SMBH model, the initial gas fraction in the disk (/gas) 
is 0.5 (0.4). The slightly lower gas content in the latter results in a 
more relaxed galaxy morphology prior to the SMBH merger, con- 
sistent with observations. 

Each SMBH particle has an initial dynamical mass of 10^* 
the total mass. The initial "seed" mass (onto which accretion is cal- 
culated smoothly via the Bondi-Hoyle-Lyttleton formula) is 7 x 10^ 
M0 for the recoil model and 10*^ Mq for the dual SMBH model. 
The BH and star formation prescriptions are described in more de- 
tail in [Springel et al.|p005) . 

In the event of a high-velocity GW recoil kick, the Bondi ac- 
cretion rate (oc Pgas/«BH) drops precipitously. Accretion onto the 
recoiling SMBH instead becomes dominated by the gas that re- 
mains bound to the SMBH. We adopt the "ejected-disk" model de- 
scribed in |Blecha et al.[p01 1^ , in which the gas ejected along with 
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Figure 1. Comparison of 9" X 9" HST/ACS image of CID-42 (a, Figure 3 from FCIO) with simulated broadband (F814W) images of the recoiling (b) and 
dual (c) SMBH models. The orientation of the NW and SE optical brightness peaks is indicated in (a); north is up, east is to the left. The simulated images 
were generated with the radiative-transfer code SUNRISE. They have the same spatial and pixel scale (0.03"/pixel) as the HST image. They are shown at the 
time of observation (tobs)> i-e. when the spatial and velocity offsets of the SMBH match the observed values. The GW recoil model shown has tijj = 1428 km 



and tr, 



- 5.9 Myr, though the galactic morphology is essentially identical in the model with tij^ = 2470 km s and 



in 



1.3 Myr. 



The SMBH is offset from the stellar density peak by 2.5 kpc. In the dual SMBH model, iobs= 



50 Myr, and the two BHs are separated by 2.5 kpc. 



the SMBH is modeled as an isolated accretion disk with a mono- 
tonically decreasing accretion rate normahzed to the Bondi rate at 
the time of the kick. 

Assuming a minimum Eddington ratio /agn ~ 
AfAGN/AfEdd for the recoiling SMBH to be observable as 
an AGN, then the observable lifetime of the SMBH is given by 



tAGK 



_3_ Mbh h 

16 MEdd /i 



BH>CJ 



Edd jEdd.mrg 



/agn 



Edd ,mrg 



-16/19 



(1) 



where ./mbh-cJ is the initial mass of the ejected disk and ,/Edd,mrg is 
the Eddington ratio at the time of SMBH merger. (A/bh and Afedd 
are time-dependent but change negligibly over the timescales with 
which we are concerned.) 

In order to compare our hydrodynamic models more directly 
with observations, we have employed the three-dimensional, poly- 
chromatic dust radiative transfer code SUNRISE j Jons son et al.| 
[2010) . The HST/ACS (F814W) image of CID-42 is shown in Fig- 
urefTt, along with SUNRISE images in the same band for our GW 
recoil model (Figure [TJ)) and dual SMBH model (Figure [TJ)). Re- 
call that in the GW recoil scenario, only one (post-merger) SMBH 
is present in CID-42, as opposed to the two SMBHs present in the 
latter model. Figures [TJ) & c show the mergers at the time of ob- 
servation (fobs), i.e. when the spatial offset (and for the recoil, the 
line-of-sight (LOS) velocity offset) of the SMBH(s) match the ob- 
served values. Although the recoiling SMBH position is marked in 
Figure [T]}, the ejected-disk luminosity is not accounted for in the 
image, so there is no observable excess of emission from this point. 
In the dual SMBH model (Figure [TJ:), the SMBH positions corre- 
spond to the two central brightness peaks. 

The galactic morphology of each model is broadly similar to 
the observed morphology. While these models are not unique solu- 
tions, they have some generic properties. For example, the promi- 
nent asymmetric tidal feature can arise in the late stages of a major 
but unequal-mass merger, such as the q — 0.5 mergers shown here. 
Fainter tidal features can also be seen in the lower-right (SW) edge 
of both the observed and simulated images. Some differences be- 
tween the two models can be seen as well. The merger remnant 
with dual SMBHs appears to be somewhat "rounder" than does the 
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Figure 2. Time evolution of relevant quantities throughout the galaxy 
merger for the recoiling SMBH (black curves) and dual SMBH (red curves) 
models. The quantity t— tmrg is plotted on the x-axis; each model has 
slightly different f mrg ■ From top to bottom, the panels show the total SMBH 
accretion rate, the global star formation rate (SFR), the total gas mass, and 
the 3-D SMBH separation prior to SMBH merger. In each panel, the dot- 
dashed red line denotes t^bs for the dual SMBH model, the black dashed 
line denotes fobs for the recoiling SMBH model with V]^= 1428 km s^^, 
and the solid blue line denotes tmrg- 
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remnant with a recoiling SMBH. However, tiie bolometric lumi- 
nosity of the recoil model (Lboi ~ 2.5 x 10^^ L©) is closer to 
that inferred for ClD-42 from SED fitting (Lboi ~ 2.9 x lO" Lg, 
FCIO). In the dual SMBH model, Lboi is only ~ 6.5 x 10^° L©. 
Table[T|summarizes the relevant physical parameters of our models 
in comparison to those inferred from observations of CID-42. 

Fig. |2] shows the merger evolution of various quantities for 
both the recoiling and dual SMBH models. Because the SFR of 
CID-42 inferred from SED fitting is ~ 25 Mq yr~^, the models re- 
quire a starburst that more or less coincides with t^hs- Similarly, the 
estimated bolometric (I/24Mm-uv) luminosity of the SE source, as- 
suming it is a Type 1 AGN, is 7.3 x 10^" L©, or ~ 4% LEdd for an 
estimatedSMBHmassof 6.5x 10^ Mq (FCIO). Thus, the SMBHs 
are also near their peak luminosity at the time of the SMBH merger. 
The presence of a stellar bulge in the progenitor galaxies acts to sta- 
bilize the gas and stellar disks against perturbations, such that the 
strongest gas inflow is delayed until the galaxies' final coalescence 
(e.g., |Mihos & Hemquist|1996| l. A high initial gas content in the 
progenitor galaxies (40-50% of the disk mass) is also required for 
fueling rapid star formation and SMBH growth. 

The dual and recoiling SMBH models have very similar evolu- 
tion, but because tahs in the dual SMBH model necessarily occurs 
some time prior to the SMBH merger (tmrg —50 Myr in our case), 
the system is observed before the peak SFR and Mbh- Specifically, 
the SFR at iobs in the dual SMBH model is ~ 6 M© yr~^, versus 
~ 21 Mq yr"'^ for the recoil model. At iobs, ~ 50 - 70% of the gas 
has been consumed. The remaining gas mass, ~ 4 - 5 x 10® M©, 
might be detectable, and an estimate of its neutral fraction obtained, 
via CO or 21cm observations. 

Our simulated merger remnants are consistent with being lu- 
minous infrared galaxies (LIRGs), as is CID-42. However, our 
SUNRISE spectra are more IR-dominated than the observed SED; 
we find that at the time of observation, -L24Mm/^boi ~ 0.2 and 0.4 
for the dual and recoiling SMBH models, respectively, compared 
with L24Mm/iboi ~ 0.1 for the Seyfert 1.8 SED fit to CID-42 
(FCIO). In the recoiling SMBH case, this disparity owes partly to 
the fact that the post-recoil AGN luminosity is calculated analyti- 
cally after the simulation, and thus does not contribute to the cal- 
culated SED. Similarly, in the dual SMBH case, the (smaller) SED 
mismatch arises partly because Lboi for the brighter AGN is a fac- 
tor of ~ 7 lower than iboi.BH inferred for CID-42 (FCIO), such 
that both AGN contribute relatively little to the calculated SED. At 
the same time, the SFR at tobs in the dual SMBH model is more 
than a factor of three lower than both the SFR in the recoil model 
and that inferred for CID-42 (EC 12). However, the results presented 
here do not depend on the shape of the simulated SED. 

2.2 Recoiling SMBH Model 

In our recoil models, the central galactic escape speed at the time 
of the SMBH merger is 773 km s~^. Thus, we expect Hk 2> v^sc 
for a recoil kick of ~ 1300 km s^^ in CID-42, such that the recoil 
event must have occurred recently. In order to match the position 
of the SE relative to the NW source, we direct the kick at an angle 
of 225° in the plane of the sky. The kick velocity then has two free 
parameters: the fraction of the total kick directed in the line of sight, 
and the amount of the measured 1360 km s^^ BL/NL offset that is 
due to the motion of the SMBH. Specifically, if the NL emission 
originates from the NW source and the BL emission from the SE, 
then a spatial-offset correction must be applied to the measured 
velocity offset, such that the actual BL/NL velocity offset would be 
1180 km s~^ (FCIO). 
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Table 1. Relevant parameters of recoiling and dual SMBH models for CID- 
42, along with the parameters inferred from observations (FCIO, FC12). Vk 
denotes the total kick velocity imparted to the merged SMBH for the recoil 
models; all other quantities are calculated at the time of observation (fobs)- 
The galactic Lboi and L24,jm for each model are obtained from radiative 
transfer calculations with SUNRISE. For the dual SMBH model, Mbh and 
Lbol.BH are given for each SMBH. The A/bh and Lboi.BH estimated for 
the SE optical source in CID-42 (FCIO) are given in the last coluirm. 

We consider models for two different kick speeds. The "low" 
kick, «k = 1428 km s~^, corresponds to (i'LOs(imrg)/i'k)^ ~ 

0. 9 and ULOs(iobs) ~ 1180 km s^^, and the "high" kick, Uk 
= 2470 km s~^, corresponds to («LOs(imrg)/'yk)^ = 1/3 and 
WLOs(tobs) = 1360 km s~^ (see Table [TJ. Note that although 
Wk ^ Wesc, there is a small difference between ULOs(imrg) and 
i'LOs(iobs) owing to deceleration. We choose these values of Uk be- 
cause velocities lower than ~ 1400 km s^^ would require the kick 
to be oriented almost entirely along the LOS, and velocities much 
higher than ~ 2500 km s^^ would have a very small probability of 
occurring. These probabilities are discussed further in § |3.1[ but we 
note that both the "low" and "high" kick velocities are quite large 
in absolute terms. 

For the low kick, tobs occurs 5.9 Myr after the recoil. The mass 
of the ejected gas disk is ~ 1.4% Mbh, and Lboi,BH = 56%LEdd 
at the time of the kick. At tobs, the AGN luminosity has declined to 
1% LEdd, or 3.2 X 10® Lq (assuming a radiative efficiency of 0.1). 
For the high kick, a gas mass of only ~ 0.5% A^bh remains bound 
to the recoiling SMBH, but because the system would be observed 
sooner (tobs = imrg +1.3 Myr), the AGN luminosity is slightly 
higher These values are summarized in Table[T] While these AGN 
luminosities are lower than that estimated for CID-42 via SED fit- 
ting (FCIO), the Eddington ratios of the order of ~ 1% are consis- 
tent. 

The relationships between tobs, tAGN, and Uk are illustrated in 
Fig. [3] For the low-luminosity AGN definition (Lboi.BH > 3 x 10® 
L0), the AGN lifetimes of the recoiling BHs (as defined in Eqn.[T|( 
are 2.3 Myr and 6.4 Myr for the low and high kicks, respectively; 

1. e., both have Iagn > tobs- In fact, for the observed parameters 
of CID-42, virtually all allowed values of i;k and fobs result in a 
low-luininosity recoiling AGN with Iagn > fobs- 

Moreover, we expect the recoiling AGN lifetimes to be even 
longer than these estimates, because the AGN luminosities and 
SMBH masses calculated self-consistently from our recoil simu- 
lations are lower than the values inferred for the SE source. We can 
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Figure 3. tobs. 'he post-recoil time at whicli tlie velocity and spatial off- 
sets match the observations, and (agn, the post-recoil AGN lifetime, are 
plotted for various total kick velocities {V]g) and LOS velocity components 
(f LOS (tmrg )/f k)- Here, fobs (solid black line) assumes no deceleration of 
the SMBH, a spatial offset of 2.46 kpc, and a line-of-sight velocity offset at 
tobs. f LOS (*obs)j of either 1360 km s~^ (top panel) or 11 80 km s~^ (bot- 
tom panel), as discussed in the text. i'LOs(*obs) is indicated by the dashed 
horizontal line in each plot. The solid black points indicate different values 
of '''LOs(*mrg)/fk> in intervals of 0.1. The colored solid lines correspond 
to t AGN for three definitions of the minimum observable AGN luminosity 
(see plot legend), using the values of A/BH(tmrg) and A/BH(imrg) cal- 
culated self-consistently from the merger simulations. Here too, values of 
''LOS (tmrg )/''k are indicated with open diamonds, in intervals of 0. 1 . The 
colored dashed lines show tAGN in a similar manner, except MBH(tmrg) 
and MBH(*mrg) are calculated such that their values at tobs match the 
observed values for CID-42. Note that while the brighter AGN phases are 
generally shorter than tobsi the low-luminosity AGN lifetimes extend be- 
yond taha for almost all vy^. Finally, the solid red triangle in each plot cor- 
responds to the kick given, where fobs assumes no deceleration. Because 
a small amount of deceleration occurs, the actual t^ha and f (tobs) are de- 
noted by the open triangles. 



use the fully-analytic ejected disk model to calculate the post-recoil 
AGN luminosities and lifetimes for arbitrary values of MBH(imrg) 
and A/BH(tmrg)- Specifically, we can match MBH(4obs) to the ob- 
served luminosity of the SE source in CID-42 and extrapolate back- 
ward to obtain its value at tmrg for each recoil model. The esti- 
mates of iAGN using these values (along with i\fBH inferred for 
the SE source) are indicated by the dashed lines in Figure|3] Here, 
the SMBH in each recoil model has I/boi,BH/iEdd > 3% at toha- 
By the lowest-luminosity AGN definition, tAGN is 4.9 Myr and 18 
Myr for the low- and high-kick models, respectively. Therefore, the 
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Figure 4. Map of gas column density from side view for our best-fit re- 
coiling SMBH models, shown at tobs- In each case, the black & white 
dot denotes the position of the BH, and the "eye" symbol denotes the ob- 
server's LOS. The color scale shows the projected gas density in units of 
cm~^. Top panel: GW recoil model with Dk = 2470 km s~^and t^bs 
= tmrg-f 1-3 Myr. Bottom panel: recoil model with tik = 1428 km s~^and 
tobs = tmrg + 5.9 Myr. 

fuel supply carried along with the ejected SMBH should not be a 
limiting factor for the recoiling AGN scenario. 

Another key consideration is whether the recoiling AGN is 
able to ionize gas and produce NL emission as it leaves the galaxy. 
Based on the current data, it is unclear whether the observed NL 
emission is coming from the SE or NW source. However, this will 
be constrained in the near future via individual HST spectra of each 
source (Civano et al. 2012c, in preparation), so it is important to 
know whether a recoiling AGN might be capable of producing the 
observed emission lines. 

It is clear from Figure |4] that the recoiling SMBH is ejected 
to a low-gas-density region within at most a few Myr. In both re- 
coil models, all of the gas near the SMBH at tobs is warm and 
diffuse, with an average temperature of a few x 10^ K and densi- 
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ties of 10~^'^ and 10~^'^ cm~^ for the high- and low-kick models, 
respectively. None of this gas is condensed into the "cold" phase 
of the multiphase model for the interstellar medium (ISM) used 
in GADGET-3 jSpringel & Hemquist|2003^ . Thus, we would not a 
priori expect any cold clouds to be locally ionized by the recoiling 
SMBH. However, sufficiently luminous AGN may ionize gas from 
outside the galaxy, similar to what is observed in some merging 
galaxies containing quasars (e.g., |da Silva et al.|201 \\ . 

[Blecha et al.| j2012[ l have developed a semi-analytic model 
for estimating the narrow-line emission around AGN in hydrody- 
namic simulations of galaxy mergers. We have applied this model 
to our simulations of recoiling SMBHs in CID-42, including accre- 
tion from the ejected gas disk. In order to estimate the lifetime of 
observable NL emission, we assume a ratio of the [O III] A5007 lu- 
minosity to the bolometric AGN luminosity of ^[oiiii/iboi.BH = 
10^^. This is the ratio of the total [OIII] luminosity to the lumi- 
nosity of the SE source estimated by FCIO. In reality, this ratio 
may be higher (such that the estimated NL lifetime would be an 
upper limit), because L[oiii] is measured from the total continuum 
rather than just the AGN continuum. However, additional uncer- 
tainty comes from the fact that I/[oiii] is measured from a slit spec- 
trum, whereas Lboi.BH is measured in a circular aperture; the di- 
rection of this effect is less clear. We therefore consider the ratio 
10~^ to be a reasonable guess, but not necessarily an upper limit. 

The NL model of |Blecha et"arij2012| l calculates the H/3 rather 
than the [O III] luminosity, because the former is less sensitive to 
the conditions in the ISM. Thus, we must also assume an [O III]/H/3 
ratio. The ratio estimated from the spectra of CID-42 is ~ 6 (JC09, 
FCIO), so we use a ratio of 10 as a conservative upper limit. This 
implies a minimum observable I/Hys/iboi.BH of 10^*. 

In this analysis of the NLR, we must account for the difference 
between the AGN luminosities calculated self-consistently in our 
simulations and the AGN luminosity of the SE source estimated 
from its F814W magnitude. Taking the values from the simulation, 
we find that the recoiling AGN can produce NL emission above 
the estimated minimum value of 3 — 5 x 10"" Lq for only about 
1 — 2 Myr after the kick, depending on Uk. In each case the NL 
lifetime is shorter than tobs. Note that any observed photoionization 
of NLRs must be ongoing; otherwise, NL AGN emission would 
essentially last for only a light crossing time ( < 10^ yr for a typical 
NLR). After this, [O III] would disappear on a timescale of ~ years 
via charge-transfer recombinations (destroying the large [O III]/H/3 
ratio characteristic of AGN), followed by H recombination on a 
timescale of < 10'^ yr. 

If instead we assume the recoiling AGN has A/BH(imrg)= 
6.5 X lO'' M0 and Lboi.BH = 7.3 x 10^° L© (the estimates for the 
SE source in CID-42 assuming it is a Type 1 AGN, FCIO), the NL 
lifetimes are longer, despite the stricter minimum H/3 luminosity 
of 7.3 x 10® L0. Specifically, the NL lifetime for the low kick is 
4.0 Myr. This is still shorter than tobs = 5.9 Myr, so the recoiling 
AGN in this model should not produce observable NL emission. 
In the high-kick model, however, the NL lifetime is 1.6 Myr, com- 
pared to tobs = 1.3 Myr. This suggests that a recoiling AGN in 
CID-42 could to produce the observed narrow lines, provided the 
total kick velocity is sufficiently large ( > 2000 km s^^). 

If the NL emission arises from stellar photoionization in the 
NW source, the recoiling AGN scenario does not require that the 
kicked SMBH produce observable NLs. However, the current spec- 
tra of CID-42 have AGN-like narrow-line flux ratios (JC09). This 
suggests that the recoil scenario requires either a large total kick 
velocity or an unusual star-forming region in the NW source. Al- 
ternatively, a recoiling AGN might be able to produce NL emission 



more locally if the SMBH happened to be ejected along a dense 
tidal stream or if the galactic line-of-sight gas distribution differed 
substantially from that of our model. 

Finally, we note that the total LOS gas column density to the 
position of the BHs is low; log A'^gas ~ 20.8 — 20.9 cm~^. Thus, 
attenuation is insignificant along the LOS to the recoiling AGN, 
even though it is recoiling away from the observer (as inferred from 
the redshifted broad H/3 line). 



2.3 Dual SMBH Model 

In the dual SMBH model, the projected SMBH separation of 2.5 
kpc occurs at fobs = fmrg —50 Myr. The slightly less dissipative 
merger in this model is key to reproducing the observed morphol- 
ogy prior to the SMBH merger. If we simply examine the pre- 
merger snapshot from the recoil model in which the SMBHs have a 
projected separation of 2.5 kpc, we find that the galaxy is substan- 
tially more disturbed, with multiple prominent tidal features and a 
more compact and asymmetric core. 

The SMBH masses at tobs are 4.5 x 10^ Mq and 2.1 x 10® 
Mq. Thus, each SMBH is more than a factor of ten smaller than 
the estimated mass of the SE source in CID-42. The less mas- 
sive SMBH is actually more luminous, with M — 13%AfEdd and 
i'boi.BH ~ 8.9 X lO'' Lq. The more massive SMBH is accreting 
at only 0.8% AfEdd, which corresponds to a bolometric luminos- 
ity Lboi.BH = 1.2 X 10^ L0. Notably, this lower-luminosity AGN 
is consistent with the upper limit on the X-ray luminosity for the 
NW source of CID-42 (Lx < 10'*^ erg s~^), so this alone does not 
place constraints on whether it may contain a SMBH. 

However, because the more luminous AGN in our model is 
about eight times fainter than the luminosity of the SE source in 
CID-42, the AGN luminosity ratio in our model may be more in- 
formative than the absolute luminosities. In CID-42, the measured 
F814W luminosity ratio of the NW and SE sources is ~ 2.5, and 
the NW source is brighter (FCIO). Given the uncertainty as to 
whether an AGN exists at all in the NW source, this tells us lit- 
tle about the possible AGN luminosity ratio. Instead, we compare 
the limit on X-ray emission from the NW source with the bolo- 
metric luminosity ratio of our simulated AGN. Luminosity ratios 
greater than 25 correspond to the 4% upper limit on the fraction 
of hard X-ray flux originating in the NW source (EC 12). In the 
kpc-scale phase of our simulation (defined loosely as SMBH sepa- 
rations < 10 kpc), the AGN luminosity ratio is above 25 for ^1% 
of the time, and not at all for separations < 5 kpc. At tobs, the lu- 
minosity ratio of the two AGN is 7.5. This suggests that an AGN 
pair with an intrinsically large luminosity ratio is a possible but rel- 
atively unlikely explanation for the lack of X-ray emission from the 
NW optical source of CID-42. 

We must also consider that a highly obscured AGN could be 
present in the NW source. Given the lack of detection of X-rays 
from this location, an obscured AGN here would have to be Comp- 
ton thick (log A'^h > 24 cm^^). The total LOS gas column densi- 
ties to the SMBHs in our model are log iVgas « 23.8 & 22. 1 cm~^ 
for the more and less luminous AGN, respectively. Notably, the 
column density of the more luminous AGN is near the Compton- 
thick limit, though these numbers should not be taken at face value. 
The Compton-thick limit refers to the equivalent neutral hydrogen 
column density of the gas, which depends on the ionic fraction of 
each element, whereas we obtain the total gas column density from 
our simulations. Our column densities therefore represent an upper 
limit on A^h- Further, the fact that the more intrinsically luminous 
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AGN is hosted in tlie more obscured density cusp works against tlie 
need for a large emergent luminosity ratio in the hard X-ray band. 

Although we do not find strong evidence for an AGN obscured 
by the galactic environment of CID-42, many AGN are thought 
to contain dense, obscuring material on small scales, well below 
the resolution limit of our simulations. If Compton-thick clouds 
fell along our line of sight to an AGN in the NW source, this 
could cause significant attenuation through the hard X-ray band. 
In this case, much of the AGN luminosity would be reprocessed 
as IR emission, creating an IR excess relative to that expected for 
a star-forming galaxy with a SFR of ~ 25 Mq yr^^. Thus, the 
fact that the SED of CID-42 is well-fit by a template for a star- 
forming galaxy (FCIO) argues against the presence of a Compton- 
thick AGN. 

Strong Fe Ka emission (EW ~ 1-2 keV |Matt et al.|2003| > is 
also expected in heavily obscured X-ray sources, (those with log 
A^H > 23 cm"^). A Fe Ka emission feature with EW ~ 100-500 
eV is indeed observed in both the Chandra and XMM-Newton data 
for CID-42 (FCIO). However, because this equivalent width could 
easily be produced by a Type 1 AGN in the SE source, we can- 
not determine whether the Fe Ko line originates from an obscured 
AGN in the NW source. 

As with the recoil scenario, we have applied the NL model of 
|Blecha et al.H2012) to our dual SMBH simulation. In this case, we 
wish to know whether the system should appear as a NL AGN and 
whether the emission lines may be double-peaked, as indicated by 
the spectra in JC09. We find that at t^ha, the lower-luminosity AGN 
has a narrow-H/3 luminosity too low to be observable. The more lu- 
minous AGN has a higher ionizing photon flux, but as it also resides 
in a higher-density region than the first AGN, its ionization param- 
eter is marginally too low by the criteria of |Blecha et al.|p012[ l. 
The lack of clear evidence for strong NL emission from these AGN 
reflects their low luminosities relative to the SE source of CID-42. 
Regardless of the NL luminosities, their relative LOS velocity in 
this snapshot is only 64 km s^^. Because this is smaller than the 
typical FWHM of NLs, it is unlikely that double-peaked emission 
lines would be resolvable in this case. Our model for CID-42 is not 
a unique solution, however, and the LOS velocities for projected 
SMBH separations of 2 - 3 kpc range from ~ — 150 km s^^. 
Accordingly, we cannot rule out the possibility that a system like 
CID-42 could produce double-peaked NLs. Civano et al. (2012c, 
in preparation) have obtained a new DEIMOS spectrum of CID-42 
with much higher S/N, which will reveal definitively whether the 
NLs are double-peaked. 



3 DISCUSSION 

3.1 GW Recoil Scenario 

Several observed features of CID-42 support the scenario that it 
contains a single, rapidly-recoiling SMBH. The offset between the 
broad and narrow emission lines is expected to be characteristic 
of recoiling AGN ( [Merritt et al.||20(j6l l. The single X-ray source 
coincident with the SE optical brightness peak is also consistent 
with the GW recoil scenario, in which the SE source corresponds 
to a recoiling AGN that has left the central stellar cusp behind 
(EC 12). Additionally, the best-fit surface brightness decomposition 
performed by FCIO consists of an extended (Sersic) component for 
the NW source and a point-like component for the SE source, con- 
sistent with a stellar cusp and a recoiling AGN, respectively. 

We note that because CID-42 is actively star-forming, imply- 
ing a high gas content, we expect the merger remnant to have a 



stellar cusp even after the displacement of stars via binary SMBH 
inspiral (e.g., |Milosavljevic & Merritt'200r;'Yu"2002','Merritt et al. 
2006} and the reco il event itself ( Boylan-Kolchin et al. 2004, Gua- 



landris & Merritt||2008^ . Gas-rich major mergers trigger a rapid 



inflow of cold gas and subsequent star formation (e.g., Barnes & 
Hernquist 1996; Mihos & Hernq uist|1996[|Mayer et al.|2007[|Hop ~ 



kins et al. 2008 1, thereby maintaining a steep central density profile. 
Stellar interactions with the SMBHs should be relatively unimpor- 
tant in such cases, because hydrodynamic processes drive the ef- 
ficient formation and inspiral of the binary SMBH i jMayer et al.| 
2007 ). Furthermore, |Blecha et al.| j20TTJ found that recoil kicks 
may actually enhance the central stellar density by displacing AGN 
feedback and delaying the quenching of star formation. We accord- 
ingly refer to a stellar "cusp" when discussing the NW source, 
though the observations do not exclude the presence of a stellar 
core on scales below the /fST image resolution of 0.03" (150 pc). 

Our recoiling SMBH models of CID-42 are able to reproduce 
the morphology and stellar mass of CID-42 in the late stages of 
a gas-rich, major galaxy merger. The peak SFR and A/bh occur 
near the time of SMBH merger, and their values at tohs are consis- 
tent with estimates from the observations (SFR ~ 25 M© yr~^ & 
iboi,BH/iEdd ~ few percent; FCIO, FC12). 

Less clear is whether a recoiling AGN could produce the ob- 
served NL emission in CID-42. If the total kick velocity is large, 
such that the system is observed soon after the kick, the recoil- 
ing AGN may still be close enough to the dense gaseous region to 
produce observable NL emission. Specifically, we find that the re- 
coiling AGN in our high-kick model (uk = 2470 km s^^) could 
produce the observed NLs in CID-42, but that it probably could not 
do so with a lower kick of 1428 km s~^. 

In the latter case, special configurations may be required for 
the recoiling AGN to produce the NL emission, such as ejection of 
the SMBH along a dense tidal stream. If the observed NLs orig- 
inate from star formation in the NW source rather than from the 
recoiling SMBH, then the recoil scenario requires unusual line flux 
ratios for a star-forming region. The origin of the NL emission is 
therefore an important consideration in evaluating the likelihood 
that CID-42 contains a recoiling AGN, and underscores the need 
for corroborating evidence. Individual spectra of the NW and SE 
sources will provide important constraints in this regard (Civano 
et al. 2012c, in preparation), and IFU data could also provide use- 
ful information about the spatial distribution of narrow (and broad) 
emission line regions in CID-42. 

IFU spectra, or HST imaging in multiple bands, would also 
constrain the colors of the NW and SE sources, and of the galaxy 
as a whole. Specifically, this would indicate whether the galaxy is 
indeed blue and actively star-forming, as predicted by SED fitting 
(FCIO, FC12) and our simulations. Such observations would also 
reveal whether the starburst is concentrated within the NW source, 
and they would provide an independent check on the starburst age 
of a few Myr estimated from SED fitting (EC 12). 

Another consideration is that a merger remnant with a recoil- 
ing AGN should contain, in principle, one central brightness peak 
(the stellar cusp) and one offset brightness peak (the AGN). An 
AGN pair, in contrast, should have both brightness peaks offset 
from the center. However, this is an inconclusive diagnostic for 
CID-42, because its center of stellar light is not well constrained. 
The surface brightness decomposition of FCIO places the center 
between the NW and SE sources, but it is consistent with being at 
the position of the NW source within the errors. The central stellar 
cusp could also be intrinsically asymmetric in the unrelaxed merger 
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remnant, and tidal features seen in projection could contribute light 
to the apparent galactic disk. 

As noted in § |2.2[ even the lowest possible recoil velocity for 
CID-42 is quite large ( > 1300 km s^'^). The intrinsic probability 
of such a large kick is not clear Recoil kick probability distribu- 
tions have been derived from numerical relativity simulations (e.g., 
[Baker et a l. 2008 , van M eter et al.|2010[|Lousto et al.|2012 1. How- 
ever, the distributions of progenitor mass ratios, spins, and spin ori- 
entations are still quite uncertain. In particular, if a circumbinary 
gas disk forms around the inspiraling SMBHs, it may align or par- 
tially align their spins (Bogdan ovic et al.|2007||Dotti et al.|2010| l. 
The maximum recoil velocity resulting from the merger of BHs 
with aligned spins is only ~ 200 km s^^. Partial alignment may 
also occur via GR precession effects l |Kesden et al.|2010{|Berti et al.| 
|2012 |, though this mechanism can also result in a«ft-alignment of 
the BH spins, yielding substantially larger kicks. At any rate, it 
is clear that the GW recoil scenario for CID-42 requires that spin 
alignment prior to BH merger is inefficient, at least in some cases. 



3.2 Gravitational Slingshot Recoil Scenario 

A rapidly-recoiling SMBH could also result from the gravitational 
slingshot mechanism. As described in § [T] this requires a long bi- 
nary SMBH inspiral time, such that a subsequent galaxy merger 
introduces a third SMBH before the initial SMBH binary has co- 
alesced. A three-body encounter may result in the ejection of the 
lightest SMBH at high speeds, and may also cause the rapid merger 
of the remaining SMBHs ^Hoffman & Loeb|[2007l l. Thus, in the 
slingshot-recoil scenario, two or three SMBHs would be present in 
CID-42: the SE source would be the ejected SMBH, and the NW 
source would contain the SMBH(s) left behind. 

While we have not constructed a model for this scenario, much 
of our analysis for the GW recoil model can be applied here. If the 
ejected SMBH retains a gas disk after the three-body encounter 
(such that it appears as an AGN), the accretion rate will evolve as 
described in § |2.1| The relationship between Uk and fobs would 
also be the same as in Figure [3] If the NW source contains an 
obscured, radio-loud AGN, this could be detected with very long 
baseline interferometry (VLBI). If it contains a binary AGN that 
survived the three-body encounter, the binary might be spatially re- 
solvable with VLBI imaging, provided its separation were at least 
a few parsec and the AGN were both radio loud. However, binary 
radio-detected AGN are intrinsically quite rare, with only one ex- 
ample known I Rodriguez et al.|2006,) despite a dedicated search of 
thousands of archival VLBI images Purke-Spolaor|2011^ . 

If SMBH binaries do not merge efficiently, the probability of 
such triple SMBH systems forming is non-negligible, especially at 
high r edshift {z > 2, jKulkami & Loeb|20I2| l. In our simulations, 
efficient merging of the SMBHs is merely an assumption owing 
to resolution limits. However, smaller-scale simulations of inspi- 
raling SMBHs indicate that gas drag (e.g., [Escala et al.pOOS) 
Dotti et al.|2007|l or non-axisymmetry of the stellar potential (e.g.. 



Berczik et al. 2006, Khan et al.|2012") can prevent SMBH binaries 



from "stalling" (the so-called "final parsec problem", e.g., |BegeLl 
[man et alT j 1980 , Yu 2002 ). CID-42 is clearly a gas-rich system with 
a disturbed stellar distribution. Thus, while the possibility of grav- 
itational slingshot recoil cannot be excluded, we consider it to be 
unlikely. Put another way, if new evidence points toward a slingshot 
recoil in CID-42 (for example, if compact radio emission were de- 
tected in the NW source and consistent velocity offsets were found 
in multiple broad lines in the SE source), the implications for bi- 
nary SMBHs would be far-reaching. This would indicate that many 



unmerged SMBH binaries may be lurking below resolution limits 
even in gas-rich environments, to the detriment of SMBH merger 
event rates for future gravitational-wave observatories. 

3.3 Dual SMBH Scenario 

The observations of CID-42 are also consistent with an SMBH pair 
rather than a recoiling SMBH. In this scenario, the NW and SE 
sources are two inspiraling SMBHs with a 2.5 kpc separation. How- 
ever, the ~ 1300 km s~^ offset of the broad H/3 line is too large to 
be associated with orbital motion on kiloparsec scales. Many AGN 
spectra show offset broad lines that are often attributed to outflows, 
but the offset line in CID-42 is atypically large. CID-42 consti- 
tutes 1/323, or 0.3%, of BL AGN in the Chandra-COSMOS survey 
with spectroscopic redshifts (El vis et al.|2009^ , and less than 1% 
of SDSS quasars have broad H/3 lines offset by > 1000 km s^^ 
( [Bonning et al.|2007[[Eracleous et al.|2012^ . 

Theoretically, a projected SMBH separation of 2.5 kpc could 
also occur during an early close passage of the merging galaxies, 
or even during a flyby encounter However, this is an unlikely ex- 
planation for CID-42. Its tidal features are indicative of a late-stage 
merger, and the highest i\'fBH and SFR typically occur near coa- 
lescence (see Fig.[2](. An unbound encounter would be required to 
explain the 1300 km s^^ BL offset via relative galactic motion, 
whereas the maximal relative speeds in our galaxy merger simula- 
tions are ~ 500 — 600 km s~^. In such a flyby, only very specific 
configurations and timing would allow a projected separation of 2.5 
kpc to be observed. 

Our dual SMBH model is well-matched to the morphology 
and stellar mass of CID-42, but a slightly less dissipative merger 
(lower /gas) is required in order to match the observed morphol- 
ogy prior to the SMBH merger, toha also occurs slightly before the 
peak SFR and SMBH accretion rate. As a consequence, the SFR, 
AIbh, and luminosities are lower in the dual SMBH model than the 
estimates for CID-42. 

Owing to these low AGN luminosities of a few x 10^ L©, their 
NL emission (as determined from the model of |Blecha et al.|20T2l > 
may be too faint to be observable. Also, their relative LOS velocity 
of only 64 km s~^ is unlikely to be spectroscopically resolvable. 
This is consistent with the supposition that the double-peaked NL 
structure in the DEIMOS spectrum (JC09) does not arise from dual 
SMBH orbital motion. 

The non-detection of X-ray emission in the NW source (FC12) 
places some constraints on whether it may host an AGN. The NW 
source could contain a quiescent SMBH with Lx < 10^^ erg s~^, 
but our models indicate that two SMBHs in a late-stage, gas-rich 
major merger are unlikely to have a large intrinsic luminosity ratio. 
Specifically, an X-ray faint AGN in the NW source requires a lumi- 
nosity ratio > 25, which occurs for < 1% of the kpc-scale phase 
in the dual SMBH simulation. 

The NW source could also contain a Compton-thick AGN. 
The upper limits on the gas column densities along the LOS to 
each SMBH (log TVgas ~ 22.1 & 23.8 cm"^) indicate that these 
AGN should be absorbed and highly obscured, respectively, but not 
Compton-thick. The hard X-ray detected AGN pair in NGC 3393, 
for example, has estimated line-o f-sight column densities of a few 
X 10^^ cm~^ | Fabbiano et al.|2()ll I. It is possible that obscuration 
on sub-resolution scales could hide an AGN in the NW source, but 
in this case, we expect the SED to be more IR-dominated than what 
is observed. Finally, while a strong iron emission line can indicate 
the presence of a heavily obscured AGN, this diagnostic is not in- 
formative for CID-42. A fairly strong Fe Kq emission line is de- 
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tected, but it is impossible to say how much, if any, of this emission 
is contributed by the NW source, which has a weaker continuum 
by a factor of at least 25. 

The detection of compact radio emission from either or both 
optical brightness peaks in CID-42 would provide strong evidence 
for the presence of one or two AGN. As mentioned above, however, 
radio-detected AGN pairs appear to be very rare, with only one 
radio-loud kpc-scale quasar pair (3C 75; Owen et^r||1985[ > and 
one pc-scale AGN pair ( [Rodriguez et al.|2006^ known. Regardless, 
given the ambiguous nature of the NW optical source and the X-ray 
non-detection there, the possibility of a radio detection is enticing, 
as this would provide the clearest evidence for an obscured AGN 
in the NW source. 



3.4 Summary 

We find that neither the recoiling nor dual SMBH scenarios for 
CID-42 can be ruled out based on current data. It is clear that 
whatever its true nature, CID-42 is not a "typical" galaxy even 
among the rather exotic subsets that host dual or recoiling AGN. 
Both scenarios require some extraordinary features, such as a dual 
AGN with an unusually large broad line offset or a rapidly-recoiling 
AGN with an atypical narrow line region. Through our analysis we 
have identified and constrained the requisite properties of CID-42 
in either case. These constraints combined with further observa- 
tions in the coming months and years, such as spatially-resolved 
and higher-quality optical spectra, radio imaging, and CO or 21 cm 
observations, should rapidly narrow the possible origins of this in- 
triguing system. 
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